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I. INTRODUCTION 



When two 2-dimensional charged systems are placed in close proximity, transport in one 
layer will drive the adjacent layer out of equilibrium. Even if the barrier separating the 
two layers is high and wide enough to prevent tunneling, interlayer interactions can still be 
sufficiently strong that a current drawn in one layer can drag along a current in the other 
layer. This phenomenon was theoretically proposed by PogrebinskiiB and by Price,i and has 
become known as frictional drag. In most frictional drag experiments a current Ji is drawn 
in one layer; the second layer is an open circuit and no current is allowed to flow. To oppose 
the dragging force, an electric field E2 develops in the second layer. The ratio of E2 and Ji 
is called the transresistivity p2i [see Eq. (|l]) below] and is a measure of a rate of momentum 
transfer from the first to the second layer. 

Experimental realizations of frictional drag between 2-dimensional systems were first 
reported by Gramila et al. for two electron layers,! and by Sivan et al. for electron-hole 
systems.0 These experiments inspired a large number of theoretical works, and the experi- 
ments (which were all done in zero magnetic field) are by now fairly well understood.i^ 

Recently, attention has turned towards frictional drag in the presence of a magnetic field 
— the topic of the present paper. Experiments of drag in a magnetic field have been reported 
by Hill et a/., El Rubel et a/.,0 Feng et a/.,lll and Eisenstein et al^ 

Frictional drag is of fundamental interest because it can serve as a sensitive probe of two 
important aspects of transport in mesoscopic systems, namely the screened interlayer inter- 
action and the form of the irreducible polarization function u;), which is central to many 
theoretical considerations. In the presence of a magnetic field, in particular, the screening of 
the interaction and the polarization function can assume different forms depending on the 
number of filled Landau levels and the degree of disorder. 

Over the past two decades new and often surprising aspects of the physics of 2- 
dimensional electron gases (2DEGs) in magnetic fields have continued to emerge. Not 
only single layer systems show intriguing physics; double-layer quantum Hall systems also 
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exhibit a number of interesting aspects,L3 and frictional drag is expected to do the same. 

In this paper we present a treatment of Frictional Coulomb drag in strong magnetic 
fields where the Landau levels are fully resolved, i.e., uJcT ^ 1 where Uc = eB/m* is the 
cyclotron frequency and r is the transport scattering time. Some numerical aspects of our 
work have been presented in a previous publication,lii here we give the full details of the 
analytic background underlying these results, extend them to other parameter values, and 
compare them critically with recent experiments.SH 

We work under the assumption that an independent electron picture applies. At suffi- 
ciently low temperatures and/or high magnetic fields this condition will not be satisfied and 
additional physics must be included as the following two examples reveal. 

In disordered systems localization becomes important at low temperatures. Chalker and 
DaniellS found that the diffusion of electrons in the lowest Landau level is anomalous, 
i.e. the 'diffusion constant' scales as D(q,uj) = Dq{uj / q^)''^^'^ , where r] is found numerically 
to be ^ 0.38. Shimshoni and Sondhi pointed out that frictional drag would be a way 
to experimentally measure 77, since the drag effect in that case would be proportional to 
T'^~^ at the lowest temperatures.0 Another example consists of high mobility systems in 
high magnetic fields, where intralayer electron-electron interactions are important at low 
temperatures. At filling factor v = 1/2 the 2DEG can be discussed in terms of composite 
fermions. The polarization function uj) assumes a unique form which was first derived 
by Halperin, Lee and Read.0 Three recent paper^Hl have considered frictional drag in this 
regime and shown that the transresistivity should be proportional to T^/^ as the temperature 
T approaches zero. 

The outline of the paper is as follows. In Section H we define the model of the system 
and establish the theoretical framework. The transresistivity is in general given in terms 
of three-body correlation functions which are considered in Section |T|. In Section ^ we 
examine in which limits the three-body correlation functions are proportional to the imag- 
inary part of the polarization function. This relation has been tacitly assumed by most 
other authors. We discuss the result and the relation to similar results in zero magnetic 
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field and a brief discussion of Hall-drag is provided. Numerical evaluations are presented in 
Section |V]; we focus on the dependence of the transresistivity upon magnetic field strength 
and temperature. Section summarizes the conclusions. 



II. FORMALISM 

We consider a system of two 2-dimensional electron gases separated by a distance d. A 
uniform, constant magnetic field B = Bz is applied perpendicular to the two layers which 
define the xy-plane. The two layers have electron densities of rii and n2, respectively. When 
a current density Ji is drawn in layer 1, the interlayer interactions will induce an electric 
field E2 in layer 2, which is an open circuit {i.e., no current is allowed to flow in layer 2). 
The induced electric field can be measured by a voltage probe, and the transresistivity tensor 
is defined according to 

E2 = P2lJl- (1) 

The transresistivity is what is measured experimentally and hence the object to determine. 
However, linear response theory, on which our theoretical approach is based, yields the 
trans conductivity (721, defined by an experiment where the first layer is biased with an 
electric field and the induced current density is measured in the second layer: J 2 = S^2iEi. 
The transresistivity can be obtained from the transconductivity by 



P21 



-0^11'?21^^22 + (712 



-P220^2lPll, (2) 



where the approximate equality is valid because we assume the magnitude of the individual 
layer conductivities, (pn)"^ and (p22)~^5 to be much larger than the transconductivity. 

To calculate the transconductivity we follow the general framework developed indepen- 
dently by Kamenev and Oreg,0 and by Flensberg, Hu, Jauho and Kinaret.ii The transcon- 
ductivity is calculated using the Kubo formula for linear response,^ i.e., it is expressed as 
a current-current correlation function 



^27(k,^^) = ^n^r(k,fi), (3) 

where Il2i'^{k, ^) is the Fourier transform of the retarded current-current correlation func- 
tion 

n^r(x-x',t-t') = -^^(t-t')([j2"(x,t) , j7(x',t')])- (4) 

Here ji and j2 are (kinematic) particle current operators in the two different layers (de- 
noted by the subscripts 1 and 2), and a and 7 are Cartesian coordinates of the two- by-two 
transconductivity tensor. (■ ■ ■) is a statistical average and [■■■,■■■] is a commutator. The 
position vectors x and x' reside in the 2-dimensional planes of the 2DEGs. In this paper, 
we assume like charges in the two layers; the sign of (J21 is reversed for unlike charges. 

To lowest order in the screened interlayer interaction, the transconductivity can be ex- 
pressed as0ii 



"^^^ ~ WA q i 2^ 



V2i{q) 



S{q,uj) 



duj 



X A2 (q, (\,uj + iri,u - z?7)A]'(-q, -q, -u - ir], -u + irj), (5) 

where is a positive infinitesimal, A is a normalization area, V21 (g) is the interlayer Coulomb 
interaction, S{q,uj) is the screening function, and hb is the Bose-Einstein distribution func- 
tion. Notice that the transconductivity tensor is a dyadic product of the two three-body 
correlation functions A2 and Ai, which are defined by 

(x, r; x', r'; x", r") = -{% {jr(x, r)A(x', r')A(x", r")}), (6) 

and will be referred to as the triangle functions. The following convention for the Fourier 
transform has been adopted 

A°(X, r; X', r'; X", r") = — ^ ^ ^ giqi-(x-x")+iq2-(x'-x")g-ia;i{T-r")~j<^2{T'-r") 

qi,q2 (^/^)^ iuii/iw-z 

xA"(qi + q2,q2,«t^i +«t^2,«tU2)- (7) 

This Fourier transform convention relies on the translational invariance of the triangle func- 
tion which applies when we consider infinite systems; extra caution should be exercised when 
considering systems of finite extent. 



III. THE TRIANGLE FUNCTION 



After the expansion in the interlayer interaction the correlation functions Ai and A2 
only depend on the individual layers and contain all microscopic details of these. To proceed 
we must choose a model for the individual layers. 

The non-interacting electron model is a good approximation for the experimental systems 
studied so faiiSiS as long as the magnetic field is not too strong and the temperature not 
too low (see Introduction). Hence, we shall treat the individual layers as non- interacting 
electrons scattering against random impurities. Within this model it has been shown that for 
short range impurity potentials, the triangle function is proportional to the imaginary 
part of the polarization function 

A(±q, ±q, ±{uj + iri), ±{uj - ir])) = q Imx(g, uj) (8) 

m* 

in absence of a magnetic field. [The same form is recovered if electron-electron scattering 
keeps the distribution function as a shifted Fermi-Diraci] Here m* is the effective electron 
mass, and the transport scattering time r is assumed to be energy independent. In the 
following two sections it will be shown that we can obtain a similar relation in the limit of 
uJcT ^ 1 for short ranged scattering potentials. 

Up to this point an ensemble average over impurity configurations has been implicitly 
understood. When the triangle function is expressed in terms of Green functions, the impu- 
rity averaging is accounted for by dressing the Green functions by self-energies and including 
vertex functions where they connect. A careful account of impurities is necessary in the pres- 
ence of a magnetic field because of the high Landau level degeneracy; leaving out impurities 
would lead to unphysical divergences. 

The free Hamiltonian is 

ff=i(P + eAf (9) 

with B = V X A. We choose to work in the Landau gauge, A = (0,5a;, 0), so that each 
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eigenstate is characterized by two quantum numbers: |A^) = |n, k). For infinite systems, the 
eigenenergies only depend on the Landau level index n\ = + l/2)huj(.. 
In terms of creation and annihilation operators, the density operator is 

p(x',r) = i E e^""-' {N\e^^-\M)ci{r)cM{T) (10) 

^ q,iV,M 

and the operator for the current is given by 

j(^) = ^ H (wiV^ + lcjv+i(r)civ(r) + W2v/ricjv_i(r)cAr(r)) , (11) 

^ N 

where we have defined the two vectors, wi = zx + y and W2 = = — ix + y with x and y 
being unit vectors defining the planes of the electron layers. When the density and current 
operators are inserted in the expression for the triangle function, we get terms involving the 
statistical average of products of 3 creation and 3 annihilation operators. The Hamiltonian 
for impurity scattering is quadratic in creation and annihilation operators, which means 
that we can use Wick's theorem to write the product of creation and annihilation operators 
in terms of products of 3 Green functions. Inserting (p!0| ) and (0) in Eq. (^) we get two 
connected diagrams. One of these is given by (see Fig. |1|) 

Zl ^ e^'i-(-'-")r(q, iV, M, r - r', r' - r")G(iV, r - r') 

^ N,N',M,q 

7"(Q = 0, N', N, t" -t,t- t')G{N', t" - r) 

r(-q, M, N\ t' - r", t" - t)G{M, t' - t") (12) 

where a current vertex 7 and charge vertices F have been included to take account of impurity 
scattering; in the other diagram the direction of the arrows is reversed. The scattering of 
electrons against impurities is evaluated in the self-consistent Born approximation. The 
self-energy diagram is shown in Fig. In the limit of lUcT ^ 1 Ando and Uemura0 have 
shown that the self-energy is given by 

he - En Fo 



he- en\ 



1 (13) 



2 2 ^ V To ; 

where Fq = (2/TT)huJc{h/T), r being the transport scattering time at zero magnetic field. 
The imaginary part of Eq. ([T3|) is taken as negative (positive) for the retarded (advanced) 



function. The width of the Landau level, 2ro, is independent of the Landau level index 
if the range of the scattering potential is smaller than the magnetic length. The choice of 
self-energy diagram implies, by a Ward identity, a specific choice of vertex functions. Born 
approximation for the self-energy implies that we must sum ladder diagrams for the vertex 
functions. In the limit of short range scattering potential the contribution from the ladder 
sum to the current vertex function can be neglected,ill i.e. 7 can be approximated by a bare 
current vertex and hence A^' = ± 1 in Eq. (|l^). For the charge vertex, on the other hand, 
the ladder sum is important. Fig. shows the diagrams corresponding to the following 
integral equation for F: 



xQ{a, tei)g{b, te2)e''"^''^'^y-'y''^^r{q, a, b, tei, le^), (14) 

where pimp is the density of impurities, Uilz) is the impurity potential, and ^ is a Matsubara 
Green function. The bare charge vertex is given by 



. -(^^/^)'yf2=(-^(^g. + g,))"-'"^r"((W2) .m<n 

fn,m{<i) = { r- (1^) 

,-(.,/2)y^2-™(-£(zg. - g,))'"-"Lr"((^g)V2) ,n<m 



where are the Laguerre polynomials and £ = yfi/eB is the magnetic length. 

In terms of dressed Matsubara Green functions and vertex functions, the expression for 
the triangle function is 

A(q, q, iQ + iuj, iu) = (wi J^i + W2 ^"2) , (16) 

where 

J-'i = F(q, n — 1, m, iui, iui + iu)T{—q, m, n, iui + iu, iui — iQ) 
X Q{n, iuoi — iQ)Q{m, iuoi + iuj)Q{n — 1, iuoi) 
+ F(q, m, n, iuoi — iu, iuJi)r{—c{, n — l,m, iuoi + iuji — iuo) 
X Q{n, iui)Q{m, iui — iuj)Q{n — 1, iuJi + iVt) (17) 
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and 



T2 = r(q, n, m, iuji, iui + 2ti;)r(— q, m, n — 1, iuji + iu, iui — iVt) 
X ^(n, iuJi)Q{m, iuJi + iuj)Q{n — 1, zcui — zfi) 
+ r(q, m, n — 1, icc;i — iu;, zco'i)r(— q, ra, m, icji + ifi, ioji — ioj) 
X ^(n, iuJi + if2)^(m, icui — iuj)Q{n — 1, icji) . 



(18) 



The summation over Matsubara frequencies iuJi can be carried out as a contour integration. 
The function A(q, q, u; + ir],u! — ir]) is then obtained by letting ifl + iu Q + u + irj, 
iQ ^ Q + if], and iu ^ uj — i77.il In the static hmit, i7 ^ 0, the result is 

A(q, q, a; + ir/, a; - ir/) = ^ ^ j ^^^(e) 



(19) 



l^wi (-P(q, rri, e, e + cj) + P(— q, n.m^e^e — uj 
-W2(P*(q, ri, m, e, e + u;) + P*(-q, rz, m, e, e - u; 



with 



P(q, rz, rri, e, e + cj) = 
r^~(q, rr, m, e, e + Ci;)r (— q, rrr, rr — 1, e + cj, e)G'^{n, e)G"'{m, e + uj)G°'{n — 1, e) 



— r (q, n, m, e, e + Ci;)r^" 
+r"'""^(q, n,m,e — u, e)T^' 
— r'^~(q, n,m,e — u, e)T^' 
+r"'""*"(q, n, m,e,e + uj)T~^'^ 
— r^^(q, n, m,e,e + u)T^' 



(-q, m, rr - 1, e + tu, e)G"(rr, e)G'"(m, e + cj)G'"(rr - 1, e) 
(— q, m,n — l,e,e — uj)G'''{n, e — uj)G^{m, e)G"'{n — l,e — u) 
(-q, m, rr - 1, e, e - uj)G'"{n, t - uj)G"'{m, t)G"'{n - 1, e - cj) 
(— q, rri, rr — 1, e + ciJ, e)G"'(rr, e)G''"(m, e + u;)G"'(rr — 1, e) 
(-q, m, rr - 1, e + u;, e)G^{n, e)G"(m, e + uj)G\n - 1, e). (20) 



The plus and minus signs attached to the vertex functions indicate the signs of the imaginary 
infinitesimals that should be added to the frequency arguments. From (|19D it can be realized 
that A is a vector with purely real components. Other general properties of A are given in 
Appendix 0. 
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IV. TRIANGLES TO BUBBLES 



We now show that in the hmit of short range scattering potentials and for uJcT ^ 1, we 
can express the triangle function in terms q and the imaginary part of the proper polarization 
function. 

The proper polarization function is obtained by analytical continuation of the density- 
density correlation function shown in Fig. ^. From the structure it is seen, that it involves two 
Green functions G, one bare vertex /, and one vertex function F; symbolically x ^ f GGT. 
The triangle function, on the other hand, involves products of three Green functions and 
two vertex functions; symbolically A ~ GGGFF. To reduce the triangle function to 
the polarization function we must therefore reduce three Green functions to two, and two 
vertex functions to one vertex function and a bare vertex. Furthermore we must introduce 
a factor of q if we want an expression similar to Eq. (H). Symbolically the task is to do the 
simplification: GG GTT — > qG G F /, which we shall now proceed to carry out. 

The key to the problem is to notice that in the expression (^) two of the Green functions 
and both the vertex functions in the product G G GTT have neighboring Landau level 
indices, n and n — 1. The retarded and advanced Green functions are given by G^'^{n, e) = 
[e - i£„ - e)]-^ Using the identity 1/AB = (l/A - l/B)/{B - A) we get 

G-(n, .)G»(n - 1, .) = t!^:;'"^°'r''f V (21) 

ujc + [n, e) — l^"-[n — 1, e) 

In the limit uJcT ^ 1 the self-energies can be neglected compared to the cyclotron frequency, 
and we can approximate 

G>..)G>-M)^^>-'>-^°'"-^-'>. (22) 

We have thus reduced products of three Green functions to products of two Green functions. 

To reduce the product of two vertex functions to one vertex function and one bare vertex 
is more involved and therefore deferred to Appendix ^ where it is shown that we can do the 
approximation in leading order of {00^7)^^ 
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F"*" (q, n, m, e, e + c<j)r {—q,m,n — l,e + uj,e)G^{n,e)G"'{m,e + u!)G°'{n — l,e) 
~ — r+"(q, n, m, e, e + u;)/m,n-i(-q)G"^(?^, e)G'"(m, e + u;) 
-— /n,m(q)r— (-q, m, n - 1, e + cu, e)G"(m, e + cu)G"(n - 1, e). (23) 

Notice that the full vertex function with same Landau indices as the Green functions is 
retained, and that the signs of the infinitesimals in the vertex function naturally follow 
the signs of the infinitesimal on the Green functions that it is multiplying. In the above 
expression there seems to be a mismatch between the Landau level indices of the Green 
functions and the Landau level indices of the bare vertices fn,m multiplying them. Matching 
indices are recovered by the identity 



Vn + l/„+i,„(q) - V^/n,m-l(q) = ^^i^/n,m(q), (24) 

which also introduces q into the expression. The square root factors in Eq. (^) come from 
the bare current vertices (see Eq. (0)). With these approximations it is a matter of simple 
manipulations to reach the following relation 

A(±q, ±q, ±{uj + ir]), ±{uj - ir])) = q x B Imx(g, uj), (25) 

which is valid for ^ L Note that for B = 0, A(±q, ±q, ±(ti; + ir]) have the same sign. 
Here the (irreducible) polarization function is 

x[G"'(m,e + t^)/„,„(q) 

X |G"'(n, e)r+"(-q, m, n, e + u, e) - G'^{n, e)r++(-q, m, n, e + u, e)| 
+G'"(m,e-cu)/„,„(q) 

x|G"'(ra, e)r (-q, n, m,e, e - cj) - G"'(r7,,e)r+"(-q, n,m, e,e - u;)| . (26) 

To obtain the transresistivity we must know the single-layer resistivities pii [see Eq. (0)]. 
For isotropic systems they have the generic structure 
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Pi: 



POi 



V 



-bi tti 



(27) 



with poi = m*/ {riie^Ti) being the resistivity in zero magnetic field. Combining Eqs. (^, 
(p5|) and (^) we find that for b ^ a the transresistivity tensor is diagonal in cartesian 
coordinates; the diagonal elements given by 



P21 



duj 



S{q,uj) 



Imxi(q,^)Imx2(q,^) 



(28) 



A semiclassical treatment of "pa yields bi = {u)cT)i, so that the term in the square bracket 
above is unity. We assume this to be the case in our numerical evaluations. As the quantum 
Hall regime is approached, bi/{uJcT)i starts deviating from 1; however this deviation does 
not change the main features of the numerical results presented in Sec. 0. 



A. A conjecture: generalization to arbitrary i3-field 

In the previous subsection it was shown, that when the individual layers are treated as 
non-interacting electrons scattering against short-range impurities, the triangle function is 
related to the imaginary part of the polarization function in the limit of uJcT ^ 1. Work has 
also been done in the small magnetic field limit uJcT <^ 1.011 We now discuss a conjecture 
for the generalization of the expression for A for arbitrary magnetic field strengths which 
extrapolates between the weak and strong field limits. 

For zero magnetic field the triangle function is proportional to the imaginary part of 
the polarization function when the impurity scattering time is independent of energy (see 
Eq. (|^)).@ One case of an energy independent transport time is when the range of the 
impurity potential is short compared to the Fermi wavelength. Likewise, for high magnetic 
fields we found that a prerequisite for a simple relation between A and Imx is short ranged 
scatterers. The task in this section is to bridge the gap between zero and high magnetic 
field. In order to do this we first observe that only two vectors can be constructed in the 
xy-plane, namely q and q x B. The triangle function is therefore of the form 
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A(±q, ±q, ±{iu + ir]), ±{uj - i7])) = A||(g, lj, B)q_ ± A±{q, B)q x B , (29) 

where the carets denote unit vectors. Knowledge of the zero and high magnetic field limits, 
results from a semiclassical analysisj^l and a perturbational calculation^ suggest the follow- 
ing conjecture for the form of the triangle function, valid for short range scatterers but for 
arbitrary magnetic field strength: 

2Titl^ ( 1 tti - \ 

Ai(±q, ±q, uj±ir],uJT iV) = 7- . , 2 ^ =F . 2 ^ ^ ^ ^"^XAQ, ^, B), (30) 

rftj y i + i + ctj J 

where ai{B) is a parameter to be determined. The magnetic field has been added as an ar- 
gument of the polarization function to emphasize that it should be evaluated in the presence 
of the magnetic field. 

The ai{B) should be chosen so that Eq. (|30D is consistent with known results. One such 
empirical result is that so far no experiment has ever observed Hall drag.ii If one assumes 
that Hall drag is absent {i.e., P21 = 0), then 

ai = — . (31) 

With the above choice one obtains in the low-field Drude limit = (co'ct)^, which is con- 
sistent with semiclassical low-field results.@ Furthermore, we note it reproduces the result 
obtained using the memory-functional formalism of Ref. 0. 

We can illustrate the plausibility of a vanishing Hall-drag by a very simple argument. 
The electrons in layer 2 are influenced by two forces which must add up to zero because 
there is no current in layer 2, i.e. 

= F12 + (-e) (E2 + (V2) X B) , (32) 

where F12 is the force from the electrons in layer 1, E2 is the induced electric field, and (V2) 
is the average velocity of the electrons in layer 2. When no current is allowed to flow in layer 
2, (V2) = 0. Furthermore, since the Coulomb force is radial, F12 must be parallel to Ji and 
hence the measured electric field is parallel to the driving current, i.e. no Hall drag. 
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However, the above argument is not as general as it seems,E3 and to second-order in 
the screened interlayer interaction, Hall-drag can occur in cases where the band-structure 
is anisotropic, which breaks in-plane inversion symmetry and when the intralayer scattering 
time is energy-dependent, which does not allow a simple description based on the polarization 
function alone. i@ Furthermore, higher-order correlation effects such as those found in an 
electron-hole system with bound pairs0 may lead to a finite Hall-drag even in the absence 
of the conditions mentioned above. 



V. NUMERICAL RESULTS 

The formula for the transresistivity Eq. (|28D must be evaluated numerically. We focus on 
the dependence on magnetic field strength and temperature. As a model for the dielectric 
function, we adopt the random phase approximations in which 

E{q, u) = [l- uj)Vu{q)f - ^)V2i{q)]\ (33) 

where Vii(g) is the intralayer- and V2i{q) is the interlayer Coulomb interaction (we have 
assumed that xi = X2 = i-e. identical layers). 

The polarization function enters both directly in the Eq. ( P5D and indirectly through 
the dielectric function. The general expression for x(g, tu) is given in Eq. (^). To make 
the numerical evaluations tractable, it is necessary to make an approximation for the vertex 
functions, which in general are given by the integral equation (|14D- In Appendix |B| we show 



that when the Landau levels are clearly resolved, we can approximate 

■n/ • I • • \ fn,m{^) /o^\ 

r q, n, m, le + lu, le) = /ofcN2T^ ^7^7 — ■ , ■ \rf ^ ^^^^ 

1 — {To/2h)^I[q, n, m)y{n, le + zuj)y[m, le) 



with T(g, n, m) given by Eq. (P^). This approximation is consistent with the assumptions 
under which we derived Eq. (^) and makes the numerical evaluation tractable. 

In order to avoid unphysical jumps in the chemical potential we must improve the self- 
consistent Born approximation (which leads to a vanishing density of states outside the 
Landau bands), and we use a Gaussian the density of states derived by Gerhardt 
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exp 



-2 



( 



£ — 



m 



) 



2' 



(35) 



To 



where a denotes spin. The chemical potential fi{B,T) is determined implicitly by requiring 
the density n to be given by 



This model has a finite range of magnetic fields where the density of extended states at 
the Fermi energy is suppressed, simulating the effect of localized states between the Landau 
bands, which are needed to obtain quantum Hall plateaus with a finite width. However, 
the quantitative details of localization, such as the critical properties of the metal-insulator 
transition, are not included in this simple model. 

With these approximations we evaluate the transresistivity given by Eq. ( pH]) as a function 
of magnetic field and temperature. For simplicity we consider two identical electron layers 
of densities ni = ?t,2 = 3 x 10^^ m~^ corresponding to a Fermi temperature of Tp ^ 120 K. 
The center-to-center distance, d, is chosen to be 800 A and the well-widths are taken to be 
200 A. The Landau level width is dependent on the transport scattering time, r, which we 
determine by choosing a mobility, er/m*, of 25 m^/Vs. The temperature dependence of the 
scattering rate - which for simplicity is neglected in what follows - will eventually lead to 
a violation of the requirement UcT ^ 1, and would set the upper temperature limit of the 
validity of the numerical evaluations. 



We focus on a magnetic field regime where the Landau levels are fully resolved. For 
simplicity we neglect spin splitting. (Because of spin degeneracy, note that filling factor u 
= odd number corresponds to a half-filled Landau level, whereas filled Landau levels have 
u = even number.) Experimentally there is a large regime where doubly occupied, clearly 
distinguishable Landau levels can be observed. For half-filled Landau levels, the density of 




(36) 



A. Magnetic field dependence 
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states g is enhanced over the 5 = value: g = gQ^2uJcT/7r where go = m* /Tth^. This is 
due to the large degeneracy of the Landau levels and implies that there are more available 
states close to the Fermi energy for the electrons to scatter into. Consequently a general 
enhancement of the transresistivity should be expected in a magnetic field. Experimentally 
the transresistivity has been found to increase as the square of the magnetic field as long 
as the Landau levels are not resolved.El When the Landau levels get resolved the picture is 
more complicated. At even filling factors the density of states is suppressed; an excitation 
gap develops, and the transresistivity should vanish as a result. These two expectations are 
both based on considerations of the density of states, i.e. the phase-space available for the 
interlayer e-e scattering. 

The screening of the double layer system is strongly affected by the density of states and 
thus also strongly dependent on the magnetic field. As the density of states at the Fermi 
level becomes smaller the electron layers lose their ability to screen and hence the effective 
interlayer interaction is enhanced. 

The resulting transresistivity can qualitatively be understood as a product of the avail- 
able phase-space and the effective interaction. In Fig. H we plot \V2i{q) / £{q,uj)\^ and 
[Imx(g, sinh(/ic<j/2A;BT)]^ as a function of filling factor together with the product of the 
two functions for a given q and uj. The maxima of the product occur at magnetic fields for 
which the filling factor is slightly above or below an even integer (where an integral number 
of Landau levels are filled). 

Fig. ^ shows the transresistivity as a function of magnetic field. At odd filling factors 
P21 is enhanced (by a factor of ~ 100 at z/ = 3 depending on the temperature) over the zero 
field value as expected. As the magnetic field is changed from an odd filling factor towards 
an even, we find that the transresistivity increases before it eventually gets suppressed when 
the chemical potential enters the excitation gap. This unusual behavior is explained by the 
competition of available phase-space and effective interaction. 

When comparing this theory with experiments, one faces the complication of spin split- 
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ting which is present in real systems. Thus, a double peak in an experimental may 
be due to two partially overlapping single-peaked structures; this is the interpretation of 
Ref. 0. However, Rubel et alU have shown experimentally that there is a regime of mag- 
netic fields (z/ = 6-15) where the single-layer longitudinal resistivity shows no spin splitting 
while the transresistivity has a clear twin-peak structure. On the other hand, their data 
at higher magnetic fields includes spin resolved structures that do not show the predicted 
double-peak structure. An improved theory, which includes spin-splitting would clearly be 
desirable. 

B. Temperature dependence 

We will discuss two regimes of temperature which show interesting behavior and which 
yield information about the polarization function and the effective interaction. 

From general properties of density-response functionsH it follows that for larger than 
the smallest relevant length scale ( = elastic mean free path for i? = 0, and I for sufficiently 
large B) and uj smaller than the inverse scattering time, the polarization function assumes 
a diffusive form 



= -g{ii)-^^ — — , (37) 



Dq^ — ioj' 

where D is the diffusion constant. For high magnetic fields the magnetic length £ = ^JhJeB 
is smaller than typical interlayer distances d {i = 180 A for B = 2 T). The dominant contri- 
bution to the g-integral in (Eq. ^8]) comes from q < 1/d, and the w-integral is dominated by 
contributions form hu ^ ksT. Hence, if the thermal energy is smaller than h/r, the diffu- 
sive form of X prevails and we should therefore expect p2i ~ In T as shown by Zheng and 
MacDonald.0 Numerical evaluations showing the In T dependence of p2i were presented 
in a previous publication.El The temperature below which the diffusive behavior prevails is 
given by ksTdm ^ h/r and is therefore sample specific. For our choices of parameters, we 
find numerically that the diffusive behavior sets in at T = 0.4 K. The T^lnT is a direct 
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consequence of the diffusive form Eq. (^^, wliicli emerges from Eq. (p6D by virtue of the 
use of the self- consistent Born approximation for the vertex correction F. As mentioned in 
the Introduction, other temperature dependences are conceivable depending on the filling 
factor, the temperature regime, and the mobility of the sample. 

For temperatures higher than T^m, the dominant contributions to the cj-integral come 
from uj > 1/t. In this regime both the real and imaginary part of the polarization function 
are strongly frequency dependent; consequently the same is true for the effective interaction 
V2i{q) /£{q,uj). Wu et al^ have studied the collective modes, i.e. zeros of Re£{q,uj). The 
absolute value of ^) falls off as a function of frequency on a scale given by the width 



of the Landau levels, Tq/H = y2co'c/7rr. For half filled Landau levels Tq/H is the frequency 
range over which the 2DEG can respond to an external perturbation. As the polarization 
decreases with u the effective interaction gets enhanced, and competition between these 
two effects leads to a non-trivial temperature dependence in the same manner that led to a 
non-trivial dependence on the magnetic field. 

In Fig. ^ we show plots of p2i/T'^ as a function of T for z/ = 3 and z/ = 5. In contrast to 
the zero magnetic field case, p2i/T'^ shows a maximum at a peak temperature Tpeak- ixhe 
peak temperature is highest for the highest magnetic field (smallest u). The maximum in 
P2i/T'^ can be associated with excitations of states where the effective interaction is strong, 
i.e. where cj) is small. As pointed out above, cj) falls off over a frequency scale 
proportional to \^ which explains why Tpeak ('^ = 3) > Tpeak(^ = 5). This prediction 
is in agreement with measurements of Rubel et a/.0 If P2i/T'^ were calculated using the 
static version of the screening function, it would be a monotonically decreasing function of 



temperature as shown in Ref. 19 



Having looked at the peak temperature for different odd filling factors, we now consider 
small changes of the filling factor around a given odd value. Specifically, we examine v = 
3 ± Su. As the filling factor moves slightly away from an odd value, the system becomes less 
susceptible to perturbations; the polarization function falls off with frequency over a smaller 
scale. As a consequence, the screening function has a minimum at smaller frequencies which 



in turn implies that one would expect the peak temperature to become smaller. In Fig. ^ we 
plot the transresistivity as a function of temperature for three different filling factors. The 
insert shows the peak temperature as a function of filling factor. We find that Tpeak indeed 
has a (broad) maximum around v = 2.8. The deviation of 0.2 away from = 3 is due to 
the general trend that Tpeak increases with magnetic field (cf. previous discussion). 

VI. CONCLUSION 

To lowest order in the screened interlayer interaction, the transconductivity of a pair 
of coupled two-dimensional electron gases is expressible in terms of three-body correlation 
functions, A, called triangle functions, which depend on the microscopic details of each 
system. In this paper we have shown that for an isotropic system of non-interacting electrons 
scattering against random, short range impurities, the triangle function is proportional to 
the imaginary part of the polarization function (Eq. p5|)) in the limit uJcT ^1. In this 
limit, we find that the transresistivity tensor is diagonal. Including bandstructure effects, 
sufficiently energy-dependent intralayer scattering time,i@ or correlations between the layers 
(in addition to the drag force) may introduce a nondiagonal elements to the transresistivity 
tensor {i.e., Hall component to the drag). 

By numerical evaluations we have illustrated how the interplay between the screened 
interlayer e-e interaction and the phase-space available for scattering leads to non-trivial 
behavior of the transresistivity as a function of both magnetic field, where the character- 
istic is the twin-peak structure, and temperature dependence, where P2i/T'^ should have a 
maximum at a temperature related to the width of the Landau levels. 

The results presented above are based on a relatively simple model for the polarization 
and screening functions. We argue that this model is applicable as long as an independent 
electron picture describes the individual layers, and should in that regime give qualitatively 
correct results when the Landau levels are fully resolved. 

Whereas the specific models for the polarization and screening functions break down at 
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higher magnetic fields and/or lower temperatures, the general expression for the transcon- 
ductivity in terms of the triangle functions remains valid (in the absence of interlayer corre- 
lations) and is open to improvements. 
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APPENDIX A: PROPERTIES OF A 

Since 0-21 = 0) must be real, A for each layer must be purely real. (To show this, 
assume that Ai has an imaginary component at qo and ujq. Then, for a purely real A2 
and |Vi2(q, = C 5(q — qo) S{uj — uq) (T21 is not purely real, leading to a contradiction.) 
Furthermore, A for each layer must be gauge invariant, since all operators which make up 
A are gauge invariant. 

Below, we give four symmetry properties of A. By definition, A(x, r; x', r'; x", r") = 
A(x, r; x", r'; x', r") which immediately implies from Eq. (|^), 

A(q, q;uj + irj,uj - ir]] B) = A(-q, -q; -u + irj, -oj - ir]; B). (Al) 

Since A is a vector quantity, in an isotropic system it must have the form A(q, q; cj + zr^, tu — 
ir]) = A||(g, uj + iri,u! — irj, i?)q + A^(g, uj + iri,uj — irj, -B)(q x B). From this, one can glean 

A(q, q;u! + iri,u! — irj; B) = — A(— q, — q; uj + iri,u! — ir]; B). (A2) 

and 

A||(q, q;LU + iri,uj - ir]; B) = A||(q, €[; u + ir] , u - ir]; -B), 

A_L(q, q; + 2?7, - ir]; B) = -Aj_(q, q;, uj + ir] , uj - ir]; -B). (A3) 
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Finally, the Onsager relationship (T21 (B) = crf;^(— B) implies 

A(q, q]uj + iri,uj — irj] B) = A(— q, — q; —u — ir], —uj + irj; — B). (A4) 

Using the above four relationships, we know how any inversion ±q, icj, ±7] and ±B affects 
A. 

APPENDIX B: CHARGE- VERTEX FUNCTIONS 

We first provide an approximation for the charge- vertex function which will also be useful 
for the purpose of later numerical evaluations. 

In the self consistent Born approximation which we have adopted for the self-energy, 
the Landau level indices are not mixed, i.e. the Green functions remain diagonal (this is 
equivalent to saying that the Landau levels are clearly resolved). Consistent with this, we 
can neglect coupling between Landau levels in the vertex functions, i.e. in the summation 
over a and h in Eq. (0) we set /fe,m(q) = /fe,m(q)5fe,m and /„,a(-q) = /n,a(-q)'5n,a- For short 
range scatterers U{k) is a constant and can be taken out of the integral. Then 

r(q, n, m, ie + iu, ie) = /n,m(q) 

-|- (ro/2/i)^X(g, n, m)Q{n, ie + iuj)Q{m, ze)r(q, n, m, ie + iu, ie) (Bl) 

with 

I{q,n,m) = (-l)"+'"e^(^'?)'/2^;^-"((£g) V2)L;;-'^((£g) 2/2). (B2) 

The vertex function is a sum of a bare vertex and a correction; we write Eq. ( PT| ) as 

r^=^(q, n, m, ei, €3) = /n,m(q) + ^''''^(q, n, m, ei, €2) (B3) 

which defines 5^^(q, n, m, ei, £2). We now show that the correction, 6, is small as compared 
to the bare vertex unless n = m. From Eq. ( |B1| ) we see that this amounts to showing that 

{To/2hfG'{n,e)G''{m,e + uj)I{q,n,m)<t:l , n ^ m. (B4) 
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We have chosen to consider F"*" (q, n, m, e, e + u) as an (important) example. For e ~ e„ = 
(n + l/2)uc and for |a;| <^ tUc we can approximate (from Eq. (pT])) 

2^ 1 

G^(n, e)G'^(m, e + < — , (B5) 

Fo [n - mjuJc 

where we have used that |G(n, e)| is of the order 2/i/Fo at its maximum. We thus have to 
verify 

To 1 



2huJr in — m 



-X(g,n,m)<l. (B6) 



In Fig. H we plot the function T(g, n, m) as a function of m for typical q and n, and conclude 
that the correction to the bare vertex is only appreciable when the Landau level indices of 
the two incoming Green functions are equal. When n = m in Eq. (|B1|) , on the other hand, 
the correction 5 is crucial for small frequencies, uj < 1/t. In this regime the correction to 
the bare vertex is responsible for the diffusive behavior which leads to a unique temperature 
dependence of the transresistivity as we discuss in Sect.|V B. 



We now proceed to explain why ( P^ ) makes the approximation in Eq. ( pB]) valid. From 
the left hand side of Eq. (El) we have terms of the form (after using Eq. (|2^)) 



F+"(q, n, m, e, e + u;)F""(-q, m, n - 1, e + cu, e)G''{n, e)G''{m, e + cu). (B7) 

Since the correction, 5, to the bare vertex is only appreciable when the Landau level indices, 
n and m, are equal, we only have to keep terms from (|B7|) with at most one correction. 
There are two terms with exactly one correction: 

/n,m(q)5""(-q, m, n - 1, e + u;, e)G'^(n, e)G'"(m, e + u) (B8) 

and 

/m,n-i(-q)5+"(q, n, m, e, e + cj)G'^(n, e)G'"(m, e + uj). (B9) 

We now argue that the first of these can be neglected compared to the second. The two 
terms should be summed over n and m (see Eq. ([ToD ). Hence, we should compare 
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Vn + l/„+i,^(q)(5~(-q, m, n, e + lu, e)G'{n + 1, e)G''{m, e + lu) (BIO) 

n,m 

and 

X! v^/m,n-i(-q)5+"(q, n, m, e, e + u;)G"'(ra, e)G''(m, e + uj), (Bll) 

where we have shifted the sum over n in ( [B10| ) in order to make the corrections, 5, directly 
comparable in order of magnitude. Since the product of Green functions is small when the 
Landau level indices are not equal, it is clear that the first term can be neglected when the 
sum over n and m is carried out. Hence the term in Eq. (|B7| ) is approximately 

r+-(q, n, m, e, e + ^)/™,„_i(-q)G"^(r2, e)G%m, e + co). (B12) 

Other terms work out similarly. 
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FIG. 1. One of the two diagrams contributing to the triangle function given by Eq. (|^. The 
Green functions (sohd hues) are indexed by quantum numbers characterizing the Landau levels 
in the Landau gauge. The propagators are 'dressed' by interactions with impurities (as shown in 
Fig. |^(a)); consistent with the Ward identity charge vertex functions F and current vertex functions 
7 are included. Note this diagram excludes some negligible contributions (see Fig. 4(a), Ref. 25). 
We assume short-range scatterers, which implies that 7 can be neglected and hence L = it 1. 
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FIG. 2. The impurity scattering is accounted for by the self-consistent Born approximation 
illustrated in (a). Thin lines are bare Green functions and thick lines are dressed Green func- 
tions. The dashed lines symbolizes interactions with impurities (crosses). Consistent with the 
self-consistent Born approximation, the charge vertex function should be taken as a ladder sum 
shown in (b). 
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FIG. 3. The density-density correlation function. Solid lines are dressed Green functions and 
the shaded area is the charge vertex function. The proper polarization function is obtained by 
doing analytical continuation, iu; — > w + irj. The explicit expression is given in Eq. ([2^). 
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FIG. 4. The two ingredients 'effective interaction' and 'pfiase-space' as a function of filling 
factor for fixed (g,a;) = (0.2/A;f, O-OOlep/^) where ef is the Fermi energy. The temperature is 
given by T/Tp = 0.01. The dashed curve is the square of the effective interaction \V2i{q) / £{q,u})\'^ 
in units of = (7r^^/m*)^, the dotted curve is 'the phase-space' [Imn(g, a;)/sinh(^a;/2fcBr)]^ 
in units of (lO^o)^- The solid line is (80 times) the product of the two. 
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FIG. 5. The transresistivity as a function of filling factor for temperature T = 1.2 K and 
interlayer distance d = 800 A. The density of the electron gases is ni = n2 = 3 x 10^^ w? and the 
mobility is 25 m^/Vs. Spin splitting has been neglected. The transresistivity shows a twin-peak 
structure: as the filling factor is changed from an odd value (where the highest Landau level is 
half filled) towards an even value, the transresistivity goes through a maximum before it gets 
suppressed. (p2i has a maximum of 1.16 J7 at = 1.8 — out of the range of the plot). 
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FIG. 6. The scaled transresistivity P2i/T'^ as a function of temperature for filling factors 
ly = 3 (solid line) and v = b (dashed line). Other parameters are as in Fig. |. /92i/r2 which m zero 
magnetic field is expected to be a constant, shows a maximum as a function of T in intermediate 
magnetic fields. The enhancement arises at a temperature which is related to the width of the 
Landau levels as explained in the main text. 
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FIG. 7. The scaled transresistivity P2i/T'^ as a function of T for three fihing factors, v = 2.6 
(dotted hne), v = 3.0 (sohd hne), v = 3.3 (dashed hne). Other parameters are as in Fig. |. The 
temperature Tpeak at which P2i/T'^ has a maximum depends on v. Tpcak as a function of v is 
plotted in the inset. There is a maximum around v = 2.8, i.e. just below an odd filling factor. 
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FIG. 8. The function I{q,n,m) as a function of m for n = 1 (circles), and n = 3 (triang 
In both cases iq = 0.35. 
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